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We study the volume averaging of inhomogeneous metrics within GR and discuss its 
shortcomings such as gauge dependence, singular behavior as a result of caustics, and 
causality violations. To remedy these shortcomings, we suggest some modifications to 
this method. As a case study we focus on the inhomogeneous structured FRW model 
based on a fiat LTB metric. The effect of averaging is then studied in terms of an 
effective backreaction fluid, ft is shown that, contrary to the claims in the literature, the 
backreaction fluid behaves like a dark matter component, instead of dark energy, having 
a density of the order of fO~^ times the matter density, and most importantly, it is gauge 
dependent. 

Keywords: Inhomogeneous Models; Averaging in Cosmology, LTB Metric. 



1. Introduction 

The amazing success of the FRW model of the universe has for years overshad- 
owed the fact that we have devised it for a smoothed out geometry of the actual 
universe which shows inhomogeneity at different scales. The dark energy concept 
cannot be properly understood until the effect of inhoniogeneities in the observa- 
tional parameters and the role of geometrical averaging is understood. While studies 
of inhomogeneous models are progressing (see Refs. [1] [2]for extensive literature re- 
view, and also Refs. [3l-|4]), the question of how to average the geometry is still an 
open question. Is it possible to write the Einstein equations for an inhomogeneous 
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universe, make then an averaging of the geometry, and obtain an effective FRW 
model of certain type? What would then be the difference between the result of 
averaging and a AC DM model of the universe ISMI? 

There have been different attempts to answer these questions. The inhomo- 
geneities maybe considered as a perturbation to the FRW mode ls; aii averaging 
process then leads to backreaction and a modified FRW universe I 2 | 7 | 8 | 9 | 10 | ll | 12 | 
The perturbative approach is likely to diverge due to the growth of perturbation at 
the epoch of structure formation, i.e. exactly the epoch of interest related to dark 
energy, as has been shown in (for a revie w of different approaches see Ref. [2|. In 
the non-perturbative approach ^^^HSIlSITT] spatial volume averaging is devised to 
smooth out the inhomogeneities of the geometry as well as the fluid content of the 
universe leading to a non-standard homogeneous FRW model. Although the meth- 
ods are non-perturbative, the difference between the real universe and the averaged 
one is usually called backreaction too. 

In this paper we are dealing with the Buchert's non-perturbative approach to the 
averaging problem in general relativity, which is based on the averaging formalism 
within the Newtonian gravity^^. While the Newtonian case, being non-relativistic, 
is well defined, the general relativistic case is to be applied with care^^. The main 
critique to the non-perturbative procedures is the fact that any inhomogeneous cos- 
mological solution leads quickly to singularities making the formalism invalid 
Using flat LTB cosmological models, we will see in detail how the singularities affect 
the volume averaging procedure. Based on resulting insights, we use the structured 
FRW model of the universe (SFRW) 20 to propose and apply two different modifi- 
cations of the volume averaging methods along the lightcone. It is then shown, that 
the averaging and the resulting backreaction is coordinate dependent, corresponds 
to dark matter instead of dark energy, and its corresponding density is 4 to 5 order 
of magnitudes less than the mean density of the universe. 

Section 2 is devoted to flat LTB metrics, a short introduction to the structured 
FRW model of the universe, and a discussion on the place of the past light cone and 
the SFRW singularities within it. After introducing the volume averaging method 
within GR in section 3, we continue with a critique of the Buchert's averaging 
method elaborating its shortcomings, suggest some modifications to it, and calculate 
the backreaction in SFRW model for different methods. We conclude then in section 
4. 

2. The structured FRW model of the universe 
2.1. Overview 

The structured FRW (SFRW) model proposed recently t^Qlis a suitable case to study 
different averaging methods in general relativistic cosmological models. The basic 
idea in developing the SFRW model is the mere fact that the events outside the past 
lightcone of the observer have had no influence on the past events observable to us. 
One may then ignore the possible inhomogeneities outside the lightcone and model 
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the universe outside the hghtcone as a FRW one. It then has to be matched to an 
inhomogeneous universe inside the past hghtcone. As the structures are effectively 
influencing the universe much later than the last scattering surface, SFRW model 
is applied to the matter dominated and pressure-less universe, and up to those z- 
values corresponding to the length scales of at least 1000 Mpc in which the universe 
is almost homogeneous. In contrast to the familiar usage of LTB metric to represent 
over- or under-density bubbles in the constant time slices of the universe, in SFRW 
the LTB junction to FRW is adapted along the past hghtcone. 



xio"= (s) Bang time (t^) 




0.1 1 10 100 1000 10000 

r (MPc) 

Fig. 1. t„ vs r for i„ = al{r^ -I- + 1). 



Although, FRW and LTB could in principle be any of the three cases of open, 
flat, or close, it has been shown in Ref. [21] that the only meaningful matching of 
these two spaces along a null boundary is the flat-flat case. Given the cosmological 
preferences, we are therefore assuming both metrics to be flat. As a consequence of 
this matching along the past hghtcone, the observed density of the universe along 
the past hghtcone has to be set equal to the density of FRW metric everywhere out- 
side the hghtcone 121!. The inhomogeneities are therefore just in our neighbourhood 
within the lightcone up to distances of the order of 1000 Mpc. 

We have to stress however that SFRW is just a toy model to understand the 
impact of nearby inhomogeneities on our observations inside our past lightcone; one 
is free to use it as a base for a Swiss cheese model, or add extra perturbations to it. 
It should not be considered as a model for structure formation similar to the onion 
model 
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2.2. The Flat LTB solution 



According to the SFRW model, a flat FRW universe outside the past lightcone is 
matched to a marginally bound or flat LTB metric inside the lightcone. Therefore, 
we restrict ourselves in this paper to the marginally bound LTB case. These are 
solutions of Einstein equations described by the metric 



R dr^ 



R'^ir,t){d9^ +sin^dd(t>'^), (1) 

in which overdot and prime (will thereafter) denote partial differentiation with 
respect to t and r, respectively. The corresponding Einstein equations turn out to 
be 

2GM{r) 



R ir,t) 
47rp(r, t) 



R 

M'{r) 
R^R' ' 



(2) 
(3) 



The density p(r, t) is in general an arbitrary function of r and t, and the integration 
time-independent function M{r) is deflned as 

M(r)=47r/ pir,t)R^dR = —p{r,t)R^ , (4) 



where p, as a function of r and t, is the average density up to the radius R{r, t). The 
metric ([T]) can also be written in a form similar to the Robertson- Walker metric. 
The following definition 

R{t: r) 



a{t,r) 



brings the metric into the form 



ds' 



'dt' 



1 + — 
a 



dr^ + r^dn^ 



(5) 



(6) 



For a homogeneous universe, a doesn't depend on r and we get the familiar 
Robertson- Walker metric. The corresponding field equations can be written in the 
following familiar form 

where we have introduced Pc{t) = 6Af (r)/r^ indicating a quasi comoving r- 
dependent density. These are very similar to the familiar Friedman equations, except 
for the r-dependence of different quantities. The solutions to the field equations can 
be written in the form 



i?(r,t) 
a(r, t) 



'9GM{r) 



■Gp,{r) 



(8) 
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These are by now standard definitions of the flat LTB metric as a solution of Emstem 
rofji 

equations 

Furthermore, the expansion and shear are defined in the following way: 



haP = gap +UaUf3 ; 

For the LTB metric (P) we obtain 



(9) 

(10) 



R' R R 
R/ 'R'R 




o- = - — - 



(11) 

(12) 

(13) 
(14) 

(15) 
(16) 



Now, the flat LTB part of the SFRW model is defined by its bang time which is 
assumed to be 



1 



(17) 



where r is scaled to 100 Mpc, as a typical inhomogeneity scale. For r-values much 
larger than 100 Mpc the bang function i„ tends to zero and, therefore, the LTB 
tends to a FRW metric. The special form of the bang function has been chosen such 
that there should be no singularity on the past lightcone. The constant a = lO^^sec 
is chosen such that the age of the universe at the last scattering surface corresponds 
to that of the standard cosmological model The comoving coordinate r is used 
alternatively in the scaled form or not, and the reader may simply see from the 
context which one is meant. Figure [T] shows i„ as a function of r. The bang time 
is almost everywhere zero except in our vicinity, reflecting the desired feature of 
SFRW. It can also be seen that for large r corresponding to the redshifts z > 1, we 
have almost FRW. Now, for this bang time we have 



0. 



(18) 



Therefore, there is no weak singularity at the origin In fact, for the LTB 
domain with this bang time, there is no singularity in the vicinity of the lightcone, 
as it is shown in figure [21 No invariant of the metric has a singularity within the 
domain of our interest. All of the quantities appearing in the averaging process 
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Fig. 2. Lightcone (time vs logio(r)) for t„ = a/(r'^ +r-^ + 1) ,a = 10^'^ sec. The solid line 
curve shows the shell focusing area R = a,t t = t„. On this curve both R' and R" are singular. 
All singularities are well inside the curve t = 3t„(r). The R' singularity (corresponding to t = 
t„ — ^'rt'„) are not shown because they almost lie completely on the curve for t = t„(r). The 
dash-dotted curve corresponds to Kretschman singularity {t = 3t„{r)). The curve plotted with 
cross signs corresponds to zero of R"; None of the geometrical scalars become singular on this 
curve. 



behave regularly at t = 3t„ where the Kretschmann invariant is singular. The 
singularities of LTB metric, which are more sophisticated than in the case of the 
Robertson- Walker metric, have been discussed extensively in the literature (see for 
example Ref. [22l - [23]) . Vanishing of each of the metric functions and its derivatives 
R, R' , R,R,R , R may lead to different singularities. In a general LTB metric there 
is another singularity, the event horizon, related to zero of 1 + where E{r) is the 
energy function of the LTB metric absent in our flat LTB case. The place of different 
singularities are summarized as follows: 

R', R", R,R' = oo — > t = tn, 

R = — > t ^ t„ , 

r'=0 t = tn + ^rt'^, (19) 

R'=0 t = tn~^K, 



'\rf f" -\- (\f f' — rf' 



2 



It is obvious that for r ^ 1, the bang time approaches a constant, and for r ^ 1 



July 23, 2009 0:15 WSPC/INSTRUCTION FILE Averaging'revised 



7 

it approaches a constant, in fact zero, meaning that for large r we have effectively 
FRW metric again. 

1.4 
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1.1 
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Fig. 3. R is plotted versus r for different time values. Looking at R' as an effective scale factor, 
it shows that the scale of the universe increases with time, although the rate of cosmic expansion 
is different at different places. 

The zero's of R are also sketched in the mentioned figures. The corresponding 
curve intersects the past lightcone on a point which corresponds to a local maximum 
of i? ' as can be seen from figure [H and none of the metric invariants, including 
that of Kretschman, have a singularity at this event. Therefore, R =0 does not 
correspond to a metric singularity. 

Having established the fact that singularities disturbing the averaging process 
happen around i = t„, we note that the time constant slices with no singularities 
corresponds to redshift values less than z = 1.45. Therefore, any averaging over 
a fixed domain includes singularities except for z-values less than 1.45. Note that 
even in domains of no singularity, where the averaged values may be defined, the 
domain gets extended beyond the lightcone and we have to face superluminal eflFects. 
Finally, we stress again that for this special choice of the bang time the condition 
(flS]) is valid and therefore, there is no weak singularity at the origin l22l. 

3. Averaging in cosmology 
3.1. Volume averaging in GR 

The cosmic fluid is assumed to be perfect and irrotational. A flow-orthogonal folia- 
tion of spacetime, i.e. synchronous slicing, with the metric ds^ = —dt^ + gijdx''dx^ 
is then assumed. We should note that the synchronous coordinates, being not nec- 
essary for the following deflnitions, are however suitable for our problem. One could 
equally choose Newtonian gauge which is more suitable to study backreaction of 
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inhomogeneities considered as perturbation to FRW models. 

The volume-averaging is based on the following definition. Let /(r, t) be an 
arbitrary function of coordinates. Its average is defined by 

= (20) 

where dV is the proper volume element of the 3-dimensional inhomogeneous domain 
D we are considering and Vd is its volume. The space- volume average of the function 
/(r, t) does not commute with its time derivative: 

(./)■- (/) = (.m-(/)w, (21) 

where the expansion scalar 9, being equal to the minus of the trace of the second 
fundamental form of the hypersurface t — const., is now a function of r and t. The 
right hand side trivially vanishes for a FRW universe because of the homogeneity. 
The averaged scale factor is then defin ed u sing the volume of our domain D by 
au = Vd (t) ^ ■ Now it can be shown that ^1^1 

9d^{0) = ^=3^= 3Hd, (22) 
V aD 

where we have used the notation gd = dao/dt, and denoted the average Hubble 
function as Hd- Therefore, all the derived quantities should be based on the average 
value aD- This is why one usually takes the above definition for the mean Hubble 
parameter and not (a/a), which is different from cld/cld- A similar difference holds 
for the second derivative of a: 

(23) 

a I (a) aD 

Therefore, the definition of the averaged deceleration parameter is not without 
ambiguity, specially because there is no firm relation like ((22)) for the deceleration 
parameter. This has motivated many authors so far to make the following definition 
for the deceleration parameter: 

aDaD aD 1 

ID = - — JTT- (24) 

a^ aD J^D 

Now, the averages of the Einstein equations, using the Ham ilto nian constraint and 
the Raychaudhuri equation, is written in the following form^ISl 

^^X='-^iPn + P.+Pn), (25) 
a_D / d 

an AttG , , , 

— = ^(Pi5+4ps), (26) 

aD 3 

where we have set {p) — pd, Pt. — S/SttG, p-n — —TZd/^Gt^G, TZ-p is the back- 
reaction term due to the average of three dimensional Ricci scalar, and S is the 
backreaction term corresponding to the non-vanishing shear defined as follows: 

E ^ {a') -I {{9- {e)r) = {a') I [{9') 91] , (27) 
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where a is the shear scalar. It may be written more suitably in the form 



where 



A 



9l] 



2,Hl = A + 3Hl, 



Substituting from (12)-(16) we obtain 



A 



R'R i?2 



R'R i?2 



(28) 



(29) 



(30) 



Note that in the synchronous gauge we have chosen, the so-called dynamical backre- 
action t errn va nishes and S, the kinematical backreaction, is equal to the total back- 
reaction H^lll, We use the backreaction here for the geometric term (G^jy) — G^^, 
where Gf^i, is the Einstein tensor and the averaging is understood to smooth out 
the inhomogeneities. Therefore, the Einstein tensor G^i^ is divided into a "homo- 
geneous" part corresponding to the FRW universe, and a remaining part which 
could be brought to the rhs of the Einstein equations as backreaction of the inho- 
mogeneities and interpreted as the energy-momentum tensor of a "geometric fluid" . 
Furthermore, the total density p, and the backreaction density pb, may be defined 
in the following way: 



P 
Pb 



PD 

Pt. - 



- Pb, 
Pn 



-1 
8^ 



(31) 
(32) 



According to this definition, the field equations lead to the following effective back- 
reaction pressure: 



1 



Pb = Pt. 



;Ptz- 



(33) 



Although 9d and Hjj are proportional, (0^) and {H'^) are not. Hence, ((28|) and ((29l) 
cannot be written in terms of H. In the case of vanishing the three-dimensional 
Ricci scalar, as it is the case in the flat LTB, the equation of state reduces to 



Ph = Pe ; = Pt,- 



(34) 



Therefore, in the flat LTB case the backreaction is defined effectively by an ideal 
fluid having the equation of state 



Pb = Pb = PT 



w — 1. 



(35) 



It is clear that it can not be interpreted as dark energy. We have not verified 
the gauge dependence of the above equation of state, as it is beyond the scope 
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of this paper. It may be different in Newtonian gauge, which would be another 
sign of coordinate dependence of the averaging method and its consequence for the 
backreaction. 



3.2. Critique of the averaging method and modifying suggestions 

The averaging process so far described and used in literature has already been 
criticized in Ref. 1191 Here we review different shortcomings of the method: 

i) It is generally gauge dependent. One then has to check each time how far the 
result in a specific gauge is viable. Note that the Eqs. (26-27) are the result of a non- 
covariant integration over the Raychaudhuri equation. This has not to be confused 
with the assumption of homogeneity of the universe and the resulting FRW universe: 
the FRW universe is a solution of the Einstein equations based on the homogeneity 
assumption. 

ii) Even if there is some natural choice of coordinates, like the comoving synchronous 
gauge, i.e. slices orthogonal to the world lines of the dust, it usually breaks down 
due to the formation of caustics 

iii) The domain of integration over which the average values are defined is fixed. In 
cosmological models of interest, even in regions without any caustics, the domain 
D extends definitely outside the lightcone for the most interesting range of time 
or redshift values, independent of how small the chosen range of r is. Integrating 
over distances outside the lightcone is, however, equivalent to taking into account 
superluminal effects. This is clearly a noncausal procedure which should not be 
implemented in the relativistic equations. It is equivalent to assigning a value to 
the density at a point on the light cone, or within it, depending on events causally 
disconnected to it, similar to the horizon problem in standard cosmology. 

3.2.1. Gauge dependence 

According to the Stewart- Walker lemma any average quantity is gauge invariant 
if the zeroth order part vanishes on the background, or is a constant scalar field there, 
or is a linear combination of products of Kronecker deltas with constant coefficients. 
The theorem has been used in Ref. [26] to show the gauge invariance of the quantities 
they are calculating. One should however be aware of the fact that the theorem is 
valid only for a perturbed metric with respect to a background. Therefore, this 
theorem can not be applied to the general case of an exact solution of Einstein 
equations such as the LTB metric we are using here. In fact, it is not clear at all 
if and when the LTB metric can be written as a perturbation to FRW metric for 
cases where backreaction maybe significant '^'^ . We will show in the section [3.2.21 
that the backreaction is non-zero for the structured FRW model in the gauge chosen, 
although it is exactly zero in another gauge indicating the gauge dependence of 
the averaging method. 



July 23, 2009 0:15 WSPC/INSTRUCTION FILE Averaging'revised 



11 

3.2.2. Modified methods: Lightcone averaging formalism 

We suggest two different alternatives to the volume averaging method in general 
relativityllS] to the aim of circumventing the shortcomings of the method elaborated 
in the section (3.2). The general motivation is to remedy the non-causal input to the 
averaging integrals through the fixed domain and also to avoid singularities. The 
first is accomplished by a time constant integral up to the lightcone. To avoid the 
singularities, we consider a region bounded by the lightcone and a time-like hyper- 
surface within it at a distance of the order of 1000 Mpc to be sure that outside this 
inhomogeneous domain we may consider the universe as homogeneous again. Note 
that the distance from any point on the curve corresponding to t = 3f„, i.e. the 
Kretschmann singular curve, is at least of the inhomogeneity order of 1000 Mpc (see 
fig [2|) . We may therefore assume that the domain within the lightcone but outside 
the Kretschmann curve is effectively homogeneous. We have already mentioned that 
the density for the events on the Kretschmann curve are regular. Therefore, these 
events are no obstacle to define the averaging if necessary. 

Both modified domains are bounded by the lightcone. To distinguish them we call 
the averaging by using the first domain the in-lightcone and the second one the 
on-lightcone averaging method. Let us summarize the three different averaging pro- 
cedures: 1) Volume averaging using fixed domain 
This is the standard averaging formalism used in literature 

2) In-lightcone averaging formalism 

The domain of integration is extended from a fixed r-value, say r = 0, to the 
corresponding r at the lightcone for a fixed time. Clearly the non-causal character 
of the standard volume averaging is absent here, but the domain may still includes 
singular points. From what has been outlined in sections 12.11 and 12. 2[ both the 
volume averaging with fixed domain and the in-lightcone averaging method may 
include caustics. Hence, these methods are only applicable to small z values where 
no singularity is seen within the domain of integration. In our model we are bounded 
to the redshifts z < 1.45. 

3) On-lightcone averaging formalism 

A time-like hypersurface within the lightcone, far from the singularities but at a 
distance from the lightcone of the order of the inhomogeneity scale such as the 
Kretschmann curve is chosen. The domain of integration is then defined for a con- 
stant time from a point on this hypersurface up to the lightcone. By this on-lightcone 
formalism we have secured the causal implications of the averaging integrals and 
also a singularity-free domain of integration. 

The first two methods are applicable along the past lightcone just up to z- values 
where the density is regular for any r value corresponding to z w 1.45. As it is seen 
from the figure [21 all the metric singularities lie under the curve corresponding to 
t = tn(r). Therefore, the hypersurface necessary for the on-lightcone method may 
be any time-like one within the lightcone bounded by the curve t = tn such as the 
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Fig. 4. Comparison of the backreaction density as a function of redshift by the three averaging 
methods for low redshifts. The value of is not defined for redshifts higher than ~ 1 due to the 
singularities in standard and In-lightcone methods. 

Kretschmann curve. The averaging integrals in this case are well-defined and maybe 
handled numerically. 

The on-lightcone averaging method is, therefore, free from both shortcomings 
and we may define all along the lightcone without incorporating any caustics or 
non-causal effects. 

3.3. Results 

Now, we have all the prerequisites to do the averaging. The results are plotted in 
the form of backreaction density and the ratio of backreaction to the matter 
density py,/ Pz versus redshift for the bang time (|17p (figures [4][6]) . Figure 4 shows 
the backreactions up to z « 1.5) for the three methods. Larger z-values leads to 
singularities and, therefore, the averaging cannot be executed in the volume aver- 
aging and in-lightcone methods. This shows again the restricted applicability of the 
volume averaging method. In the case of on-light-cone method we have been able to 
integrate up to the surface of the last scattering corresponding to z k. 1100. There- 
fore, any meaningful statement about the effect of backreaction can only be made 
for this method. For the sake of completeness, and as a reference to our numerical 
calculation, we have also calculated numerically the case of bang time t„ = equiv- 
alent to FRW with vanishing backreaction plotted in ( figures [71 [8]) . The vanishing 
of backreaction in this case is easily seen up to the numerical errors. In all the cases 
under consideration, the backreaction density and pressure are positive and roughly 
5 order of magnitudes smaller than the dark energy needed to explain the accel- 
eration of the universe. Hence, it cannot be concluded that the backreaction has 
any effect towards explaining the mysterious dark energy. The backreaction for a 
LTB universe using the volume averaging for a fixed domain (figure |4|) shows clearly 
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z 

Fig. 5. Backreaction density as a function of redshift by the on-lightcone averaging. 




1 10 100 1000 



Fig. 6. Ratio of the backreaction density to the matter density by the on-Ughtcone averaging. 

that the effect is even three order of magnitudes smaller than all the other values 
for the SFRW model. A noticeable effect is the increase of the backreaction with 
redshift after its minimum at about z « 4. This is in contrast to the motivation 
of SFRW in choosing the bang time (fT7| . One would expect that the backreaction 
would decrease all along the past lightcone. To understand this behavior, we have 
plotted the bang time (fT7|) versus the redshift in figure 9. It shows clearly that the 
bang time has a minimum at roughly the same redshift value of the minimum of the 
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Fig. 7. Backreaction density as a function of redshift by the on-lightcone averaging for t„ = 0. 

backreaction. This redshift value corresponds to r w 5000 Mpc which is of the order 
of magnitude of the inhomogeneity of the universe. It is, therefore, clear why the 
backreaction increases after that minimum: the bang time, and as a consequence 
the deviation of LTB from FRW, increases. The selected bang time (fT7|) seems to be 
well-behaved up to the inhomogeneity scale. We know already that the universe is 
homogeneous at that scale to a very good approximation, and is modeled by FRW 
standard cosmology. Hence, we may assume that within the SFRW philosophy the 
metric after z « 4 is given by FRW. 

After publishing the first version of this work, a study SI appeared, discussing an 
averaging approach along the past lightcone, in contrast to time constant domains. 
The authors define a new procedure for averaging in cosmology not related to the 
volume averaging and its modifications we have propos ed i n this paper. A separate 
study on the backreaction in a flat LTB model universe was published, in which 
the authors suggest a so-called 'running averaging scale' to modify the volume av- 
eraging method, using a gauge similar to that of Ref. [TT] and, as in Ref. [TTl show 
the vanishing of the backreaction in it. They are mainly interested in the modifi- 
cation of the luminosity distance. A recent paper on the cosmological backreaction 
from perturbations'^ also appeared. Working in the Newtonian gauge, the authors 
announce a backreaction of the order of 10~^ times the matter density, similar to 
our result, but with an effective equation of state w « —1/19. 

4. Conclusion 

We have explored the possibility of explaining, at least partially, the dark energy 
using the volume averaging of Einstein equations for a specific model based on the 
flat LTB inhomogeneous solutions. We were able to define two modifications to the 
familiar method with a fixed domain of integration without the usual shortcomings 
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Fig. 8. Ratio of the backreaction density to the matter density by the on-lightcone averaging for 

tn = 0. 




Fig. 9. Ratio of the bang time to the time on the lightcone versus redshift. 



such as the caustics and non-causal implications. It turned out that the backre- 
action in the gauge chosen mimics a normal matter with positive pressure, and a 
density roughly 4 to 5 order of magnitudes less than the matter density. This is in 
contrast to vanishing of backreaction in other gauges, and confirm the claims that 
backreaction is gauge dependent. We therefore conclude that, although the method 
can be made free of singularities and superluminal effects, the volume averaging in 
the SFRW toy model leads to a non-vanishing addition to the normal matter five 
order of magnitudes less than the matter density, indicating a kind of dark matter 
component, but no sign of any dark energy component. 
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